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ABSTRACT 



In risk assessment, it is most important to know the nature of all available information, data or model parameters. 
More often, it is seen that available information is interpreted in probabilistic sense because probability theory is a very 
strong and well established mathematical tool to deal with aleatory uncertainty. However, not all available information, 
data or model parameters are affected by aleatory uncertainty and can be handled by traditional probability theory. 
Imprecision may occur due to scarce or incomplete information or data, measurement error or data obtain from expert 
judgment or subjective interpretation of available data or information. Thus, model parameters, data may be affected by 
epistemic uncertainty. Fuzzy set theory or possibility theory can be explored to handle this type of uncertainty (Dutta et al. 
2012). Sometimes, it is also seen that some model parameters are affected by aleatory uncertainty and some other model 
parameters are affected by epistemic uncertainty then there is the need for joint propagation of uncertainties. In this paper, 
an attempt has been made to combine probability distributions, generalized fuzzy numbers and normal fuzzy numbers 
within the same framework. 

KEYWORDS: Aleatory & Epistemic Uncertainty, Fuzzy Set, Generalized Fuzzy Number, Hybrid Method, Risk 
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The aspect of uncertainty is an important and integral to any risk assessment process. For any decision making 
process involving risk the modelling and quantification of the uncertainties is required. The uncertainties are basically two 
types viz., aleatory and epistemic. Aleatory uncertainty arises due to inherent variability, natural stochasticity, 
environmental or structural variation across space or through time, manufacturing or genetic heterogeneity among 
components or individuals, and Varity of others sources of randomness. On the other hand epistemic uncertainty arises due 
to the insufficient knowledge about the world, which includes small sample sizes, detection limits, imperfections in 
scientific understanding etc. When some parameters are affected by aleatory uncertainty and other by epistemic 
uncertainty, how far computation of the risk, then one can either transform all the uncertainties to one type or use some 
methods for propagation of the uncertainties. Many researchers have studied the issue of combining probabilistic and 
possibilistic representation of aleatory and epistemic uncertainty respectively within the same computation of risk. For 
example, Guyonnet et al, (2003), kentel et al (2004) has proposed hybrid method for join handling of probability and 
possibility distributions. The hybrid method proposed in Guyonnet et al, (2003) combines the random sampling of 
probability distribution functions (PDFs) with fuzzy interval analysis on the a-cuts. In order to compare random fuzzy set 
to a tolerance threshold Guyonnet, et al (2003) performed a post-processing of this result. Baudrit et al (2004) laid bare a 
shortcoming of this post-processing method. Baudrit et al (2004) showed how the theory of evidence, also called theory of 
Dempster-Shafer (or theory of belief functions; Shafer, 1976) could provide a simple and rigorous answer to the problem 
of summarizing the results of the hybrid computation for comparison with a tolerance threshold. In the hybrid approach 
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proposed by Kentel et al (2004) combined utilization of fuzzy and random variables produces membership functions of risk 
to individuals at different fractiles of risk as well as probability distributions of risk for various alpha-cut levels of the 
membership function. Li et al. (2007) has proposed different hybrid method for joint handling of probability and possibility 
distributions. Dutta et al (2012) proposed a hybrid method to deal with both variability and uncertainty within the same 
framework of computation of risk. 

In this paper, an attempt has been made to combine probability distributions, generalized fuzzy numbers and 
normal fuzzy numbers within the same framework. 

BASIC CONCEPT OF FUZZY SETS 

Environmental/human health risk assessment is an important aid in any decision-making process in order to 
minimize the effects of human activities on the environment. Unfortunately, usually environmental data tends to be vague 
and imprecise, so uncertainty is associated with any study related with these kinds of data. Fuzzy set theory provides a way 
to characterize the imprecisely defined variables, define relationships between variables based on expert human knowledge 
and use them to compute results. In this section, some necessary backgrounds and notions (Dutta et al., 201 la) of fuzzy set 
theory that will be required in the sequel are reviewed. 

Definition 

Let X be a universal set. Then the fuzzy subset A of X is defined by its membership function ju A : X [0, 1] 

Which assign a real number ju A (x) in the interval [0, 1], to each element x e A, where the value of ju A (x) at x 
shows the grade of membership of x in A. 

Definition 

Given a fuzzy set A in X and any real number a e [0, 1]. Then the a -cut or a -level or cut worthy set of A, 
denoted by A is the crisp set 

a A = {xeX: MA (x)>a} 

The strong a -cut, denoted by a + A is the crisp set 
a A = jxe X :ju A (x)>a} 

(A) 

Definition 

The support of a fuzzy set A defined on X is a crisp set defined as 
Supp(A) = {xeX :jU A (x)>o} 
Definition 

The height of a fuzzy set A, denoted by h(A) is the largest membership grade obtain by any element in the set and 
it is denoted as 

h(A) = sup /u A {x) 

xeX (4) 
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Definition 

A fuzzy number is a convex normalized fuzzy set of the real line R whose membership function is piecewise 
continuous. 

Definition 

A triangular fuzzy number A can be defined as a triplet [a, b, c]. Its membership function is defined as: 



, a < x < b 



b — a 



,b < x < c 



Ic-b ( 5 ) 
Definition 

A trapezoidal fuzzy number A can be expressed as [a, b, c, d] and its membership fuzzy number is defined as: 



,a < x < b 

b — a 

1, b < x < c 
d — x 



d — c 



,c < x < d 



(6) 



GENERALIZED FUZZY NUMBERS (GFN) 

The membership function of GFN A= [a,b,c,d;w] where a<£><c<<i,0<w<lis defined as (Chen, 11985, 

1999) 

O, x < a 

x — a 

w ,a < x < b 

b — a 

w, b < x < c 



M A (x) 



w ,c < x < a 

d — c 

0,x > d 



(7) 



If w = 1, then GFN A is a normal trapezoidal fuzzy number A = [a, b, c, d]. If a = b and c = d, then A is a crisp 
interval .If b = c then A is a generalized triangular fuzzy number. Ifa = b = c = d and w = 1 then A is a real number. 
Compared to normal fuzzy number the GFN can deal with uncertain information in a more flexible manner because of the 
parameter w that represent the degree of confidence of opinions of decision maker's. 

PROPOSED HYBRID APPROACH 



Consider the model 

M = g (P„ P 2 ,..,P m Qj, Q 2 ,-,Qr, Fj, F 2 ,...,F n ) 



(8) 



which is a function of parameters. Suppose P h P 2 . . . P m are m parameters presented by probabilistic distributions 
while Q b Q 2 ,...,Q r are r parameters presented by generalized fuzzy numbers with heights w b w 2 , ...,w r respectively and F h 
F 2 . . ., F„ are n parameters presented by normal fuzzy numbers. 



The approach is explained below: 
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• Consider all r generalized fuzzy numbers QnQ,2,-»,Qr & n normal fuzzy numbers F h F 2 ,..,F n . Then calculate 
a -cut for each fuzzy number. As height of each generalized fuzzy numbers are different so to perform 
computation we consider a values lie within the interval [0, w], where w = min (w b w 2 ,..., w,). 

First calculate a (=0)-cut for each fuzzy number. Then r+n numbers of closed intervals (as a -cut gives closed 
intervals) i.e., 2{r+n) numbers of values will be obtained. 

• Generate m number of uniformly distributed random numbers from [0, 1] and perform Monte Carlo simulation to 
obtain m number of random numbers by sampling probability distribution. 

• Assign all m random numbers and 2(r+n) values in the model M and calculate M/ - Inf(M) and Mf up = Sup (M). 

• Repeat stepl to step3 for 5000 times 

• Plot cumulative distribution function of (M/" f , M 2 inf , ... , M 500 o" f ) and (M/" p , M 2 S " P , M 500 o" p )- 

• Calculate a -cuts for other a values (note that greatest a value is w. i.e., a e [0, w] ) of each fuzzy number and 
repeat step 2 to step 5. If proceed in this way a family of cdfs will be obtained 

From theses family of cdfs membership functions (generalized fuzzy numbers) at different fractiles can be 
generated with each of height w. 

A CASE STUDY 

To demonstrate and make use of the proposed hybrid method a hypothetical case study for cancer risk assessment 
is presented here. Suppose water became contaminated due to the release of radionuclide to the water. Need to calculate 
cancer risk for the ingestion pathway. 

The risk assessment model due to the ingestion of radionuclides in water as provided by EPA, 2001 is follows 

C xIRxEFxED 

Ra*= xCSF 

EWxAT 

(9) 

Where C is concentration (mg/L), IR is the ingestion rate (L/day), EF is the exposure frequency (days/year), ED is 
exposure duration (years), BW is the body weight (kg), AT is averaging time (equal to 70 years x 365 days/year), and CSF 
is the cancer slope or potency factor associated with ingestion (mg/kg-day) -1 . 

Values of the parameters for the calculation of cancer risk are given in the table 1 . 

Table 1: Parameter Values Used in the Risk Assessment 



Parameter 


Units 


Type of 
Variable 


Value/Distribution 


Concentration (C) 


mg/L 


Probabilistic 


Normal(0.15, 
0.0005) 


Intake rate(IR) 


L/day 


GFN 


[4.75, 5, 5.25; 0.8] 


Exposure frequency (EF) 


Days/year 


Constant 


350 


Exposure Duration (ED) 


Years 


Constant 


30 


Average Time (AT) 


Days 


Constant 


25550 


Body Weight (BW) 


Kg 


Fuzzy 


70 


Cancer slope factor 
(CSF) 


(mg/kg-day)" 1 


Constant 


[0.14, 0.15,0.16] 
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The result of the cancer risk assessment due to the ingestion of radionuclides in water using equation (1) is 
depicted in figure (1). For simple and clear representation here we consider a = 0.0, 0.4, 0.8. 



jt=0{L) 

jt=0[U) 

jt-0.4(L) 

jt=0.4[U) 

jt=0.B[l) 

- - 



Figure 1: Cumulative Distribution Functions of Cancer Risk When a = 0.0, 0.4, 0.8 

In this study, representation of the parameter concentration (C) is considered as probabilistic distribution having 
with mean 0.15 and standard deviation 0.0005, intake rate (IR) is taken as generalized fuzzy number while representation 
of the parameter cancer slope factor is considered as normal triangular fuzzy number and other parameters are considered 
as constant. Using our proposed method to deal with both type uncertainty natures in the risk assessment we have the result 
in the form of Cdfs. From these cdfs, risk at different fractiles (Maxwell et al, 1998, Kentel et al, 2004 & Dutta et al, 
2011b) can be calculated and which are obtained in the form of generalized fuzzy number. For instance, at 90 th fractile, 
cancer risk value lies in the generalized fuzzy number [5.881e-05, 6.546e-05, 6.721e-05, 7.428e-05; 0.8]. Similarly, at 80* 
fractile, risk values lie in the generalized fuzzy number [5.872e-05, 6.534e-05, 6.711e-05, 7.417e-05; 0.8]. Their graphical 
representations are given below. 



T 



Figure 2: Cumulative Distribution Functions of Cancer Risk When a = 0.0, 0.4, 0.8 
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Figure 3: Membership Function of Risk at 90 Fractile 
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Figure 4: Membership Function of Risk at 80 Fractile 
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CONCLUSIONS 

Risk assessment is an important and popular aid in the decision making process. The aim of risk assessment is to 
estimate the severity and likelihood of harm to human health from exposure to a substance or activity that under plausible 
circumstances can cause to human health (Kentel et al. 2004). In risk assessment, it is most important to know the nature of 
all available information, data or model parameters. More often, it is seen that available information model parameters, 
data are tainted with aleatory and epistemic uncertainty or both type of uncertainty. When some model parameters are 
affected by aleatory uncertainty and other some parameters are affected by epistemic uncertainty, how far computation of 
the risk is concern, one can either transform all the uncertainties to one type of format or need for joint propagation of 
uncertainties. In this paper, an effort has been made to combine probability distributions, normal fuzzy numbers and 
generalized fuzzy numbers within the same framework. To demonstrate and make use of the proposed hybrid method a 
hypothetical case study for cancer risk assessment is carried out here. After performing risk assessment using our approach 
risk is obtained in the form of Cdfs and from which, membership functions of the risk are generated at different fractiles. 
The membership functions of risk at different fractiles are generalized fuzzy number since representation of at least one 
parameter is taken as generalized fuzzy number. Also we have observed that shape of the generalized fuzzy number is 
trapezoidal type it is because any arithmetic operation of generalized fuzzy numbers (also generalized fuzzy number and 
normal fuzzy number) produces trapezoidal type generalized fuzzy number. 
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